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Abstract 

We investigate the pseudo specific heat of SU(2) gauge theory near the crossover 
point on 4 4 to 16 4 lattices. Several different methods are used to determine the 
specific heat. The curious finite size dependence of the peak maximum is explained 
from the interplay of the crossover phenomenon with the deconfinement transition 
occurring due to the finite extension of the lattice. We find, that for lattices of size 
8 4 and larger the crossover peak is independent of lattice size at f3 co = 2.23(2) and 
has a peak height of CV, co = 1.685(10). We conclude therefore that the crossover 
peak is not the result of an ordinary phase transition. Further, the contributions to 
Cy from different plaquette correlations are calculated. We find, that at the peak 
and far outside the peak the ratio of contributions from orthogonal and parallel 
plaquette correlations is different. To estimate the finite temperature influence on 
symmetric lattices even far off the deconfinement transition point we calculate the 
modulus of the lattice average of the Polyakov loop on these lattices and compare 
it to predictions from a random walk model. 



1 Introduction 



The pseudo specific heat Cy of SU{2) gauge theory was already investigated in 
the beginning of Monte Carlo lattice studies. It is known to have a peak near 
(3 = 4:/g 2 fa 2.2, in the crossover region between strong and weak coupling behaviour. 
A first finite size analysis by Brower et al. Jl| on 4 4 — 10 4 lattices revealed a strange 
dependence of the peak on the volume V = (iV CT a) 4 of the lattice. Here, a is the 
lattice spacing and N a the number of points in each direction. The location of the 
peak shifts with increasing volume from smaller (3— values to larger and then to 
smaller ones again; the peak maximum decreases with increasing volume. Such a 
behaviour is unknown for any ordinary phase transition. The nature and origin of 
the peak remained therefore unclear, though a connection to the nearby endpoint of 
the first order critical line in the (3— adjoint)-plane was proposed by Bhanot and 
Creutz 0: "The peak in the specific heat of the SU(2) model is a shadow of this 
nearby singularity". Recently, new interest in this extended SU(2) model came up 
again || [|] and a line of second order deconfinement transition points was found, 
connecting the normal deconfinement transition point at (3a = to the endpoint of 
the first order critical line inside the ((3, /3A)-plane. The origin of the change from 
second to first order behaviour and the possibility of a bulk transition have still to 
be examined. 

For both the relation to the extended SU(2) model and the unexplained finite size 
behaviour in the crossover region a new study of the pseudo specific heat is useful. 
Moreover, we may now determine Cy with much higher statistics and also on larger 
lattices than in the early calculations, and we can apply modern analysis techniques. 
In addition, since symmetric lattices are often used to simulate zero temperature 
physics, it is important to estimate remaining finite temperature effects which may 
show up in the pseudo specific heat . 
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2 Methods to Calculate C v 



We use the standard Wilson action for SU (2) 

S = (3 ■ £ P^x) , (2.1) 



X,fJ,U 



where 



P^(x) = 1 - ^TrU^x) , (2.2) 



is the plaquette or energy and U^ u (x) is the plaquette link operator. The sum extends 
over all independent forward plaquettes. There are N P = QN* such plaquettes. We 
denote the lattice average of the plaquettes by P 

P =TrT. • (2-3) 

P X,\XV 

The speudo specific heat is then defined by 

C - d{P) - B" d{P) (2±\ 
There are three methods to determine Cy : 

i) one measures the plaquette expectation values (P) as a function of f3 and calculates 
the numerical derivative at (3m — P + A/3/2 from 

Cv(Pm) = -§^((P)((3 + A/3) - (P>(/3)) ; (2.5) 

ii) one measures the variance of the plaquettes, which is proportional to Cy 

C v = P 2 N P ((P 2 ) - (P> 2 ) ; (2.6) 



or, 
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Figure 1: The pseudo specific heat Cy calculated from the variance (filled sym- 
bols), the numerical derivative (open symbols) and the DSM interpolation (lines), 
calculated on a 4 4 (circles) and a 6 4 lattice (squares). 

iii) one calculates the sum of plaquette-plaquette correlations 



C v = (3 2 £ ((P^(x)P^(x')) - (P)'< 



(2.7) 



The most straightforward way is, of course, to calculate the variance of (P). The 
density of states method (DSM) || may then be used to interpolate between the 
points. With increasing lattice size the DSM requires however more and more sim- 
ulation points to obtain a densely populated action histogram. This comes about, 
because the variance of P is essentially proportional to Np 1 ( implying a nearly 
Np— independent specific heat Cy, except for the smallest lattices ). From 



S = pNpP 



(2- 



1/2 

it follows then, that the width of an action histogram for a fixed p— value is ~ Np . 
Compared to that the size of the complete multiple points histogram for the action is 



3 



of order Np. Thus more simulation points are needed for larger lattices to interpolate 
in the same (3— range. Consequently we have only applied the DSM to the 4 4 and 
6 4 lattices. 

In Fig. m we show the results from methods i) and ii) and the DSM interpolation 
for these lattices. At each (3— value we took on the average 90-120 thousand mea- 
surements. Between the measurements five updates, consisting of one heatbath and 
two overrelaxation steps were performed, so that the autocorrelation time was of 
the order of one. As can be seen from Fig. [l] there is complete consistence between 
the different methods. It is perhaps appropriate at this point to note that if, as 
it is often done, the plaquettes are measured during and not after each update the 
resulting plaquette variance is about 20% smaller than expected, because of local 
correlations among the plaquettes, though the plaquette average is correct. 
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Figure 2: Contributions to Cy on an 8 4 lattice from parallel (circles) and orthogonal 
(triangles) plaquette correlations up to distance R as a function of R. The sum is 
plotted as squares. The results are shown for (3 = 2.23 (a) and (3 = 2.70 (b). 

We have also investigated the plaquette correlations. We find in general a rapid 
fall with R — x' — x, the correlation length is of order 1-2 lattice spacings. The pla- 
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quettes P fJ , l ,(x) and P^ u '(x') in Eq. ( |2.7| ) may be in parallel or orthogonal planes. At 
the peak ((3 w 2.23) we find that the total contribution of the orthogonal correlations 
is about 30% higher than that of the parallel correlations, whereas far away from 
the peak, at (3 = 2.70, the contributions are essentially equal. This is demonstrated 
in Fig. 0, where we compare the contributions to Cy in Eq. ( |2.7p . Shown are the 
corresponding sums up to distance R as a function of R. Note, that the maximal 
nontrivial (with respect to the periodic boundary conditions) diagonal distance R 
on an lattice is N a a. We see that at the crossover the orthogonal correlations 
reach their total contribution only at R « 4a, whereas at f3 — 2.70 a distance of 
R fa 2a is sufficient. In contrast to that the parallel correlations have a much shorter 
range and distances of R > 1.5a play no role, both at the crossover and at higher 
/3— values. There is also no difference among those parallel correlations where the 
two plaquettes are in the same plane or in different parallel planes. 



3 Finite Size Dependence of Cy and Finite Tem- 
perature Effects 

In Fig. |( a) we compare the results for Cy from lattices with N a = 4, 6, 8, 12 and 
16. The general behaviour already found in [0J is fully confirmed. On the other 
hand, we see that there is no further finite size dependence in the peak region, if 
N a > 8. This is shown in Fig. |3](b) in more detail for the crossover region. The 
results suggest, that the finite size dependence of the smaller lattices is related to a 
different phenomenon. Indeed, the critical point for the N T = 4 finite temperature 
deconfinement transition is at (3 C = 2.30, very close to the crossover peak positon. 
Since we are using periodic boundary conditions for all directions, the approach 
to the critical point corresponding to N a will influence the plaquette expectation 
values. To check this, we have calculated the Polyakov loop 

1 N a 

L(f) = -Trn^,t , (3.1) 
A t=i 



and its lattice average L 
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Figure 3: The pseudo specific heat Cy vs. /3 on lattices. Part (a) shows the 
results for N a = 4 and 6 from the DSM interpolation, for N a = 8, 12 and 16 the 
measured points were connected by lines to guide the eye. Part (b) shows the 
crossover region for N a = 8, 12 and 16. 



As can be seen from Fig. [|, the expectation value of the modulus of L is not 
zero on symmetric lattices, not even in the strong coupling limit, i.e. we have finite 
temperature effects also on symmetric lattices. Well below (3 c {N a ) the quantity 
is a constant. With increasing f3 it starts to increase already before the transition 
point. It is obvious, that due to the nearby transition points the crossover peaks 
of the 4 4 and 6 4 lattices are more distorted than those of the larger lattices, where 
only the right shoulders of the peaks are slightly influenced. 

We can actually calculate the ^—dependence of at (3 = from a simple 
random walk model. For (3 = the Polyakov loops L(x) at the N% positions x are 
a set of equal random variables. The expectation value of the modulus of their sum 
is then assumed to behave as 



6 



(3.3) 



for large N a , so that 



(\L\)p=o = c-N^ 2 , (3.4) 
where c is a constant. Indeed, a simulation at (3 = confirms this relation in detail. 
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Figure 4: The expectation value of the modulus of L for N a = 4, 6, 8, 12 and 16 vs. (3. 
The broken vertical lines show the locations of the corresponding finite temperature 
phase transitions. 

The logarithms of the (\L\)— values, which we calculated at (3 — are shown in Fig. 
|5| together with the corresponding results for (3 = 2. At this (3— value all lattices 
apart from the 4 4 lattice have already reached their respective strong coupling limits. 
Deviations from the limit indicate then the onset of finite temperature effects. One 
may thus test, in which (3— region results from a specific symmetric lattice can be 
used to simulate zero temperature physics. For SU{2) a fit to the (3 = data leads 
to c = 0.400(1) . The relation Eq. ( |3.4| ) holds as well in SU(3) gauge theory and we 
find there the value c = 0.296(2). A first estimate of the proportionality constant 
c is obtained from the expectation value of the modulus of a single L(x), i. e. the 
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result for N a — 1. For SU(2) it is 

(|L(f)|) /3=0 = 4/37r , (3.5) 
which differs from c only by 5%. 
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Figure 5: The logarithm of (\L\) vs. lnN a on 4 4 - 16 4 lattices at (3 = (solid line), 
at P — 2 ( circles ) and from an estimate from the (\L\) value at N a — 1 (dashed 
line) . 



Our final conclusion is, that the crossover peak is not the result of an ordinary 
phase transition. For large lattices the peak is at (3 co = 2.23(2), its height is C v , co = 
1.685(10). For small lattices (N a < 6) however, the interplay of the crossover 
phenomenon and finite temperature effects shift and distort the peak considerably. 
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